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We propose a novel approach to define and measure the statistics of work, internal energy and dissipated heat 
in a driven quantum system. In our framework the presence of a physical detector arises naturally and work and 
its statistics can be investigated in the most general case. In particular, we show that the quantum coherence of 
the initial state can lead to measurable effects on the moments of the work done on the system. At the same 
time, we recover the known results if the initial state is a statistical mixture of energy eigenstates. Our method 
can also be applied to measure the dissipated heat in an open quantum system. By sequentially coupling the 
system to a detector, we can track the energy dissipated in the environment while accessing only the system 
degrees of freedom. 


I. INTRODUCTION 

One of the cornerstones of the Copenhagen interpretation 
of quantum mechanics is the measurement postulate: after a 
projective measurement, the wave function collapses into an 
eigenstate of the measured observable. In this framework, two 
subsequent measurements of an observable are not indepen¬ 
dent, as the first measurement perturbs the state of the system 
and thereby affects the result of the second [1], Still, there 
are quantities in classical physics which are not ’’local” in 
time and need two (or more) observations to be determined. 
Among them are the charge flowing through and the work 
done on a system. In such cases, the extension of classical 
definitions and protocols to the quantum realm is not straight¬ 
forward. Recently, the statistics of the work done on a quan¬ 
tum system and, more generally, its energy exchanges have 
attracted much attention [2-9], Besides a fundamental inter¬ 
est, the thermodynamics of quantum systems has important 
implications to the energetic performance of quantum devices 
[10] and quantum heat engines [11]. 

An established protocol to measure work involves a dou¬ 
ble projective measurement of the energy of the system at 
the beginning and at the end of the evolution. Such a two- 
measurement protocol (TMP) can be described in terms of 
classical conditional probabilities [2, 3, 12, 13]. It has proven 
successful in formulating quantum fluctuation relations in a 
setting where the system is initially in a statistical mixture of 
energy eigenstates. However, TMP has limitations that have 
so far failed to receive the due attention. These limitations be¬ 
come apparent when one tries to apply TMP to a more general 
class of processes, namely, those in which the system is ini¬ 
tially in a quantum-coherent superposition of different energy 
eigenstates. Most quantum gates developed in the context 
of quantum information and computation [10] belong to this 
class. The problem with TMP is that the initial measurement 
forces the system into an eigenstate of the initial Hamiltonian. 
The ultimate result of this operation is to reduce the dynamics 
to a classical statistical one [7] and to destroy the interference 
effects that generate, in Feynman’s words, the ’’interfering al¬ 
ternatives” in the dynamics [14]. In this respect, TMP fails 
to answer in a general way the most straightforward and im¬ 
portant question [15, 16]: how much energy is it needed in 


order to perform a given quantum operation on an arbitrarily 
prepared quantum system? 

In this Article, we address the key question we have just 
posed by proposing a measurement protocol that is meant 
to preserve the quantum-mechanical nature of the work per¬ 
formed on a quantum system. In this protocol, a quantum 
detector is coupled to the system at the beginning and at the 
end of the evolution. The information on the energy is stored 
in a phase shift that can be measured, for example, by inter¬ 
ferometric means. Our protocol can be formally derived from 
a path-integral description of the dynamics by adding a con¬ 
straint on the admissible paths [17-19]. Its predictions coin¬ 
cide with that of TMP for a mixed initial state. However, as 
soon as we introduce quantum coherence in the initial state, 
we find a stark disagreement between the two protocols, even 
at the level of the first moment of the distribution, that is, the 
average work. We discuss the reasons for this disagreement 
and set the stage for further investigation. We also discuss 
how to extend our protocol to measure the dissipated work in 
driven, open quantum systems, by accessing the degrees of 
freedom of the system only. As compared to previous propos¬ 
als relying on measurements of the environment [6, 20, 21], 
our protocol may provide an experimentally more accessible 
way to measure the statistics of heat and work in this case. 


II. GENERAL FORMALISM 

It is known that the work done on a quantum system cannot 
be associated to a hermitian operator and, therefore, it is not 
an observable [3, 22]. In general, the work performed on an 
open system depends on the full evolution of the system and 
not only on its initial and final state. As quantum trajectories 
(or paths) play a key role in determining energy exchanges 
and dissipation, we find it natural to tackle the work measure¬ 
ment problem by using a path integral approach [14], The 
formalism we describe in this section is an adaptation of that 
developed by Sokolovski in a series of papers [17-19]. We 
refer to the Appendix A and to the original papers for more 
technical details. 

We consider a closed quantum system whose dynamics 
is generated by a time-dependent Hamiltonian Hy(t). The 
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drive starts at t = 0 and ends at t = T. The corresponding 
unitary evolution operator can be approximated as U (T) = 
exp (—i [q dtHs(t)) « where T* denotes the 

time-ordering product and in the second writing we have dis¬ 
cretized the time in A + 1 steps of length At and used the no¬ 
tation Hs(kAt ) = 6g. 

The probability amplitude to go from the initial state | i//o) 
to a given final state | yz-j) can be decomposed into Feynman 
paths [14]. Differently from the usual approach, in which the 
dynamics is described in the position-momentum basis, we 
exploit the freedom to choose any complete basis at each time 
to decompose the paths. If we are interested in the behavior of 
a time dependent operator A(t), the preferred basis is the one 
composed of its eigenstates, i.e., A(t) \cii(t)) = a,-(f) |a,-(f)). 
The idea behind this choice is that it allows us to associate 
to A{t) a value a(t) depending on the path traversed during 
the evolution. In a more formal way, we write 4 = kAt and 
{|a,-(4))} = {|a*)} for a complete basis set. By inserting 
the completeness relation for |a)) into the expression for the 
probability amplitude to go from | VA)) to | we obtain (see 
A) 

(y 7 \U(7)\y 0 ) » (ra|nf =0 e“' A ^|v/o) = £ Ay (1) 
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where CP is the path defined by the sequence of states 

(see Fig - and is the P robabilit y 

amplitude to go from i//o) to | associated to that path. 
Along the path CP, the operator A(t) takes the set of values 

,a?, = a(t). Thus, we can also associate to CP any 

functional F[CP] of a{t). 

At this point, we add a constraint to the evolution by re¬ 
quiring F[CP] to take the value /. The constrained probabil¬ 
ity amplitude reads A[f\ = J]y5(F[CP] — f)Ay. As in [19], 
we consider functionals of the form F[CP] = dtfi(t)a(t) = 
A? Lit where /3(f) is an arbitrary function. The Dirac 

delta in the expression for A[f\ can be written as a Fourier 
transform in a conjugate space described by the variable A, 
as follows: 5(F[(P] — /) = / dX exp [/A(F[CP] — /)]. Notice 
[17] that A and / can be thought of as eigenvalues of conju¬ 
gate operators A and / acting on an additional Hilbert space, 
their corresponding eigenstates |A) and |/) satisfying the re¬ 
lation (A|/) = e~A/. Denoting with A x =A(kAt) and recall¬ 
ing that Af/\a^ ) = af N \af N ), we can write (see Appendix A) 

A\f \ = f dX Ly Ayf, where 

A X J =(v*r,A |.... 

■(4| e -^F^ 0 )|a« ) . ( aO| lKO) / ) . ^ 

is the probability amplitude to go from the state |y/o,/) to the 
state | l/r T ,A) [17-19]. The evolution described by A X J is gen¬ 
erated by the effective Hamiltonian 

H(t)=H s (t)-X p(t)A(t). (3) 

Equation (3) plays a central role in our work and it is worth a 
few comments, (i) The additional Hilbert space we introduced 



FIG. 1. Quantum work and path integral. Pictorial representation 
of the unitary evolution of a quantum system from the initial state l/fo 
(in this case an eigenstate fo the initial Hamiltonian) to the generic 
final state \jf<j, described in terms of paths in energy space. The time- 
dependent energy spectrum e, (f) of the system Hamiltonian is plotted 
in black. Quantum trajectories (blue, red) consists of a sequence of 
jumps between different eigenstates. The red trajectory satisfies the 
constraint (in this case At/ = 0) while the blue ones do not. 


can be related to a detector in the von Neumann measurement 
scheme [23]. Therefore, requiring that the functional F as¬ 
sumes the value / along the evolution is equivalent to intro¬ 
ducing a detector and coupling it to the observable we wish 
to measure. Here, A and / act as the momentum and posi¬ 
tion operator of the detector, respectively, (ii) The interaction 
described by (3) does not induce any transition between the 
eigenstates |A) of the detector momentum, (iii) The informa¬ 
tion about the system-detector interaction - and hence about 
the value taken by the functional F - is encoded in the phase 
accumulated between the eigenstates |A) and |A'). 

Observation (iii) suggests that the statistics of the integrated 
observable A(t) can be determined by measuring the phase of 
the detector, as done in the full-counting-statistics approach 
(FCS) [24—26]. Let the composite system be initially de¬ 
scribed by the factorized density operator po = p® <g>pj), where 
p!) and p{) are the density operators of the system and the de¬ 
tector, respectively. Then the phase difference acquired be¬ 
tween the eigenstates \X/2) and \—X/2) of the detector reads 




(A/2MQI— A/2) 
<A/2|pO|—A/2) 


= Tr s 


Ui / 2 (t)p 0 sU f _ x/ 2 (t) 


(4) 


where Ux(t) = it exp [—/ f^dt 1 (Hs — X[5A)\ is the evolution 
operator generated by (3). The function G, plays the role of a 
moment generating function, as the «-th moment of A is given 

by (A") = (-i) n d n Gx/dX n \ x=0 [27, 28], 


III. INTERNAL ENERGY OF A CLOSED SYSTEM 

We now have the instruments to determine the variation of 
the internal energy of a driven closed system. Starting from 
Eq. (3), we take A (t) = H$(t) and /3 (/) = 8(t — T) — 8(t) [19]. 
This corresponds to coupling the detector and the system only 
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at the beginning and at the end of the drive. (More precisely, 
we couple the system and the detector at time t = 0 ~ and at 
time t = T + , i.e., immediately before and after the starting and 
ending drive.) The resulting G x is given by Eq. (4), with the 
evolution operator (see Appendix A) 

U x/2 {7) = . (5) 

The so-obtained G x is a measurable quantity and can be 
used to determine all moments of the internal-energy varia¬ 
tion A U in the system. However, the interpretation of this 
result presents some subtleties, which we are now going to 
discuss. 

It is known from previous work [26-29] that, in general, the 
Fourier transform of the G% in Eq. (4) cannot be associated to 
a probability distribution. A similar problem is encountered 
when defining the FCS of electron transfer across a supercon¬ 
ducting device [30]. If a probability distribution cannot be 
defined for the variation of the internal energy, the question 
then arises what is the meaning of the moments generated by 
G a . To clarify this point, let us first analyze the first moment, 
which for a closed system corresponds to the average work 
performed on the system. A physical expectation for the re¬ 
sult can be developed by considering the following gedanken- 
experiment. We repeatedly prepare the system in the same ini¬ 
tial state ps(0). Half of times we just measure Hg(0) and de¬ 
termine its average (fis{ 0)). The remaining times we first ap¬ 
ply the desired evolution to arrive at ps(t) = U(t)ps(0)U^(t) 
and then measure Hs(7) to determine (Hs( T)). According to 
this procedure, we estimate variation of the internal energy as 
AU = (Hs(7)) — ( Hs{0 )). This result is the same as obtained 
from Eq. (4); by contrast, it cannot be reproduced by TMP. 
To pinpoint the differences between the two methods, let us 
explicitly write AU as obtained from Eq. (4) : 



e °)+ E PhoJo^iUl, 

GoA/o 


( 6 ) 

where p°,. Jo = (e io \p° s \e jo ), U k , = (ef|f/(T)| £ ,°), u[ k = 

(e?|l/ + (T)|ef), W Ki = |<% t |[/(T)|£P)| 2 , and |e°> and \ef) are 
the eigenstates of the Hamiltonian at the beginning and at the 
end of the evolution, respectively. The first term in ( 6 ) is the 
same as in TMP [2, 3, 31] and can be straightforwardly inter¬ 
preted in terms of classical conditional probabilities. On the 
contrary, the remaining terms, which depend on the initial co¬ 
herences Pg jo ■ , are of a purely quantum nature. These terms 
are destroyed by the initial measurement of H 5 ( 0 ) performed 
in TMP. The fact that the interfering terms can have important 
effect in the statistics of the work was first pointed out in Ref. 
[15, 32], 

The situation is well exemplified by the cyclic evolution of 
a coherent superposition of energy eigenstates into itself. As 
both the initial and final state and the initial and final Hamil¬ 
tonians are the same, we would naturally expect AU = 0. 
However, this needs not be the case in TMP, as the individ¬ 
ual energy eigenstates after the first measurement can evolve 
into different states with different energies. As specific ex¬ 
ample we consider the case of two level system driven by a 


periodic Hamiltonian so that H${ T) = //$(0). We initialize 
the system in a state \jfo) that, a part from a phase factor, 
in left unchanged by the evolution generated by Hs(7), i.e., 
U (T) i//q) = ey/o). The existence of such a state is guar¬ 
anteed, for instance, by Floquet theorem [33]. Clearly, the 
internal energy of the system does not change and, therefore, 
AU = 0. This is correctly predicted by the first moment cal¬ 
culated from G). 

But, in general, | v Po) needs not be an eigenstate of Hg( 0). 
We consider the case in which | v Po) = cos a |ej) + sin a £ 2 ) 
where | £,:) (i = 1 , 2 ) are eigenstate of the initial (and final) 
Hamiltonian. If we take a to be a free parameter, then the 
requirement of cyclic evolution for | v Pq) forces the evolution 
operator to take the form in the {|£i), (£ 2 )} 


/ cos§-Mcos2asin§ isin 2 asin^ 

' l i sin 2 a sin B, cost!; — f cos 2 a sin ^ 

With the TMP, after the first measurement, the system is found 
in |£ 1 ) with probability cos 2 a and in (£ 2 ) with probability 
sitr a. These two states now evolve independently as the “in¬ 
terfering alternatives” have been destroyed by the projective 
measurement. The final result for the work distribution can 
be computed in terms of classical conditional probabilities I) y 
for the system to make a transition between states i and In 
particular, for the average change in the internal energy, one 
finds 


AU = AE(P \2 — P 21 ) = A£’cos2a sin 2 2a sin 2 2E, , 

where A E = (£ 2 |f/s( 0 )|£ 2 ) — (£i|fls(0)|£i). We have thus 
found that AU is generally nonzero, except in the cases 5 = 0 
(trivial evolution), a = 0,7 t/2 (the initial state is an energy 
eigenstate), and a = n/A (equal superposition of the two 
eigenstates). 

The interpretation of higher-order moments is not trivial 
and stands as an open question in the field [28], As already 
mentioned, it is known from the FCS [26-29] that the Fourier 
transform of G x is a quasi-probability which can assume neg¬ 
ative values. A probability distribution can be retrieved in 
some cases after partial integration of the relevant Wigner 
function [26, 27]. Ultimately, these complications are rooted 
in the full quantum treatment of the detector [27]. Indeed, 
different types of measurements performed at the end of the 
evolution yield different distributions for the same quantity, 
each of which must be interpreted accordingly. The mea¬ 
surement of the phase of the detector has the advantage that, 
since A is a constant of motion, preserves the “quantumness” 
of the evolution and leads to Eq. (4). It is our belief that the 
quantum correlations stemming from Eq. (4) should not be ig¬ 
nored; instead, they deserve further exploration. For instance, 
the negativity in the quasi-probability distribution of work can 
be thought of as due to nonclassical temporal correlations of 
the energy, leading to the violation of a Leggett-Garg-type in¬ 
equality [27, 28, 34], Further progress in this direction will 
hopefully appear in future work. 
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IV. OPEN SYSTEM AND HEAT STATISTICS 


We now turn our attention to the more general case in which 
the system is coupled to an environment during the drive. In 
order to determine the work performed on the system, we need 
to complement the measurement of the internal energy dis¬ 
cussed above with one of the dissipated heat. To this end, 
different approaches have been proposed, including the mea¬ 
surement of an engineered environment [3, 20, 21, 35, 36]. 
Yet measuring the environment is a challenging task, restrict¬ 
ing the applicability of these proposal to specific physical re¬ 
alizations. In the following, we describe an extension of our 
measurement protocol that allows one to obtain the statistics 
of the work and dissipated heat by accessing only the system 
degrees of freedom. 

We describe the open system by the Hamiltonian H = 
Hs + Hse + He where He and //<,/. are the environment and 
system-environment coupling Hamiltonians, respectively, and 
assume weak coupling between the system and environment. 
We first take both H and //<, to be time independent and con¬ 
sider a measurement of H$. Then (5) simplifies into U x /2 = 

e -i\H Se -i7H e i\H s anc j jj s 

are constant, no external 
work is done on the system and the variation of internal energy 
must correspond to the dissipated heat. We can also show (see 
Appendix B) that the statistics obtained from the above equa¬ 
tion is the same as the one obtained by measuring the environ¬ 
ment degrees of freedom [2, 35, 37, 38]. We conclude that for 
an open system with constant Hamiltonian, the scheme gives 
the statistics of the dissipated heat Q (Appendix B). 

For a time-dependent Hs(t ), we discretize the evolution 
in N time intervals At, denote H k = H l s + Hse + He with 

U k = e Within each time interval At, the Hamiltonian 

is constant. At the beginning and at the end of each interval, 
we instantaneously couple our detector to H|. In analogy with 
(5), the evolution operator for each interval reads 


Ui /2 = e 


-i&8£ -iAtft K J±H‘, 


Ae 


Xftk 


e l ^ n s 


(7) 


Each is defined so that we keep track of the heat Qt dis¬ 
sipated in the time interval ( k — 1 )At <t< kAt. As a result, 
the information on the dissipated heat along the evolution is 
stored in the phase of the detector. Notice the opposite sign 
in the exponents with respect to Eq. (3) takes into account the 
fact that an emission (absorption) by the environment, i.e., de¬ 
creasing (increasing) of the environment energy, corresponds 
to an absorption (emission) process of the system, i.e., in¬ 
creasing (decreasing) of the system energy. 

In order to account for the variation of the internal energy as 
well, we must add a measurement of Hs at the beginning and 
end of the evolution (Appendix C). Putting things together, 
the total evolution operator reads 


Ux/: 


= e‘2 ft sn"U k > 


4=o ji e 


• A tj 0 
~ l 2 M S 


( 8 ) 


A pictorial representation of the scheme described by (8) is 
presented in Fig. 2. In the case of unitary evolution, H k = 
H$ and we immediately recover the closed-system result for 



FIG. 2. Measuring work and dissipation in open quantum sys¬ 
tems. Schematic representation of the sequence of driven evolutions 
and interactions with the detector with the open system protocol in 
Eq. (8). The evolution steps exp (~iAtH k ) are represented by the flat 
line and characterized by the Hamiltonian H k . Each coupling with 
the detector is represented by a dot. The coupling is either of the 
form exp (—iXH k ) (red dots), or exp ( iXH |) (blue dots). In the blue- 
shadowed region, the evolution is frozen (Ap = 0) and the Hamilto¬ 
nian changes by AH. In the red-shadowed region, the Hamiltonian is 
constant (AH = 0) while the density operator changes by Ap. 


the variation of the internal energy. The moment generating 
function is the same as in (4) with Ux given by (8). Let us 
calculate its first moment, which gives the average work W = 
—idGx/dX |i = q. We find (Appendix C) 


W = Tr$ 


Hs Ps,N 


HsPs,o — ^kHs^-Ps,k > 


(9) 


where A p$, k = Ps.k Ps.k-\- In the first two terms we 
recognize the variation of the internal energy of the sys¬ 


tem: AU = Ti\ 


H^ Ps,n ~ H® psfi 


Accordingly, we iden¬ 


tify the remaining term with the dissipated heat: Q = ^. k Qk = 


£<.-Tr s 


H k s Ap s , k 


. Notice that while AU depends only on the 


initial and final state of the system, Q is determined by the full 
dissipative evolution as in the classical counterpart (Appendix 
C). 

In the fast-decoherence limit, the dissipated heat takes an 
illuminating form. When energy-relaxation processes are 
much faster than the dynamics induced by the drive, we al¬ 
ways find the system in its instantaneous thermal equilibrium 
state Ps(t) = exp [— Hs(t) / kgT] / Zs(t) where Z$(t) is the par¬ 
tition function of the system and T is the temperature of the 
environment. In other words, the system evolves through 
states of quasi-equilibrium. Defining the Von Neumann en¬ 
tropy as S = Tr[pslogpy], we can show (see Appendix D) 
that the variation of the density operator is related to it by 


Tr s 



ksTAS k where AS k is the variation of entropy 


at time t k . Then we can link the variation of entropy to the 
dissipated heat by the relation Q k = k/fl A S k , confirming the 
above interpretation of Q as the dissipated heat. 

There is an alternative way to interpret Eq. (9). Taking 
the time derivative of the average internal energy, we have 


d(H s (t))/dt = (Hs(t)p(t)) + (Hs(t)p(t)} [15], If the evolu¬ 


tion is unitary, the second contribution vanishes and we can 
relate the variation of the system Hamiltonian to the instan¬ 
taneous work done on the system. By expanding the product 
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in (8), we identify pairs of sequential system-detector interac¬ 
tions of the form exp (/A//£ +1 /2) exp (—fA//|/2). Each such 
pair effectively keeps track of a variation in the Hamiltonian. 
As the variation is instantaneous, the system has no dynamics. 
We can interpret the action of the pairs as a “measurement” 
of the work done on the system by an external force. This 
interpretation is strengthened by the analysis of Eq. (9). By 


regrouping the terms, we can write it as W = Try 


T. k AH^p sk 


where A= fl| +1 


-m. 


One may wonder whether the repeated coupling to the de¬ 
tector can ’freeze’ the dynamics of the system (dynamic Zeno 
effect). This turns out not to be the case: a dynamic Zeno ef¬ 
fect would require A —» °°, while we derive our physical quan¬ 
tities, i.e., the moments of the work done, in the opposite limit 
A —> 0. Our protocol can instead be regarded as a noninvasive 
measurement [28, 39] of the work distribution. In fact, the 
moments generated by G; depend on evolution operators that 
describe the dynamics of the open system without a detector. 


V. CONCLUSIONS AND OUTLOOK 

In summary, we have shown that the statistics of work per¬ 
formed on a quantum system exhibits nonclassical correla¬ 
tions in a deeper and more fundamental way that it had so far 
been appreciated. Such correlations become apparent once 
one replaces the customary double projective measurement 
with a less-invasive coupling to a quantum detector. The re¬ 
sulting protocol is immediately applicable to the case of uni¬ 
tary evolution and can be suitably extended to treat open quan¬ 
tum systems. Our approach puts the problem of work under 
a new perspective and leads the way toward further investiga¬ 
tions. In particular, the links between quantum-mechanical 
work, Leggett-Garg-type inequalities [27, 34], weak mea¬ 
surements [39], and stochastic quantum trajectories [40-42], 
await to be fully elucidated. An experimental test of our pre¬ 
dictions is in reach of state-of-the art quantum technology. 
Among different architectures, superconducting quantum cir¬ 
cuits in combination with nearly-quantum-limited parametric 
amplifiers are a first choice, given the high degree on control 
achieved in recent experiments [41, 43-45], 
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Appendix A: Probability amplitude in the path integral 
representation 


We consider a closed quantum system whose dynamics is 
generated by a time-dependent Hamiltonian Hs(t). The drive 
starts at f = 0 and ends at t = T. The corresponding unitary 
evolution operator is 


U(7) = 7V'Ac UtMt) ~ n % =0 e~ iAtA s (Al) 

where 7 denotes the time-ordering product and in the second 
writing we have discretized the time in N + 1 steps of length 
Af and used the notation Hs(kAt) = fl|. Our goal is to write 
the probability amplitude to go from | y/o) to | y/<j) in terms of 
Feynman paths [14]. Equation (Al) is approximated to the 
order Af 2 and it is a convenient way to describe the evolution 
in terms of path integral. As a final step, we will take the 
limit Af —► 0 to recover the usual continuous description of 
the evolution. 

At any time t k we can find a basis {| ai(t k ))} = {|af )} such 

that the completeness relation Y.i k 4) (d\ = 1 holds. In this 
notation, k is a time index and i k denotes the eigenstate ba¬ 
sis index at time t k . By inserting the completeness relation 
for |a*) into the expression for the probability amplitude we 
obtain 


(wmvwo) = 

k =0 4 

{a 2 h \e- lMH s | al) (a\ { \e~ ,AtH ° |a°) (a° | yr 0 ) . (A2) 

The term (a k ik \e~ l ^s |a* - j) is the probability amplitude to go 
from |a* _ J) to |af ). Then, 


at = (Yr\e- iAt "" 14:;}... 

(af \e~ ,AtH s |4)<4 \e~ ,AtH s\a\) (a-J V'o) (A3) 

is the probability amplitude to go from | yf$) to | yr<j) passing 
through the sequence of states: |a? ), |a--),..., \af }. This se¬ 
quence define a path 7 in the basis space {Ja} }}. We interpret 
(A2) as the sum over all the possible paths of the probability 
amplitudes: 


(yr T \U(7)\yr 0 )= £ A v . (A4) 

all 9 
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In the limit Af —► 0, we obtain a continuous path a(t). In this 
way, we can associate to the path 7 a physical quantity F[(P] 
depending on it. Mathematically, F is then a functional of 7. 

We now add a constraint and select only the paths that sat¬ 
isfy some properties. We are asking which is the probability 
amplitude A[f\ to go from | VA)) to |VAr) though a path 7 de¬ 
termined by a(t) for which the functional F\7\ assumes the 
values /. The constrained probability amplitude reads 

A[f\=Y J 8(F[7}-f)Ay (A5) 

where the delta function restricts the admissible paths to those 
for which F[7\ = f. F[fP] could a generic functional of 7. 
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However, as in Ref. [19], we assume that the functional de¬ 
pends on the integral of the path a(t) (or \aj t ), |a? ),| af N ) 
in the discretized expression) 

r 7 N 

F[?}= dtp ( t)a(t ) = Af £ (A6) 

■ /0 t=o 


where [3(f) is an arbitrary function. 

By writing the Dirac delta in terms of Fourier Transform 
5(F/P] —/) = f dX exp [— iX{f — F[CP])] inserting it in the 
path integral representation, and splitting the term F [T J ] with 
respect to the corresponding time interval, we obtain 


A[f] = / dX 


e~ iXf Z (VT\e- iAtA se imi}N <\a? N )... 


/ 1 | i&hPoaV i 0 w 0 i„, \ 

\ a h\ e Se °\ a i 0 ) \ a i 0 \VO) ■ 


(A7) 


We must choose the time-dependent basis set {\a\ }} consid¬ 
ering the observable A(t) we are interested in. In particular, 
we must take them in order that A & \a\) = cjA \a\ ). For small 
Af [46], we can write 


A[f \« / dX 

J rp 


(al\ e - iAt ^-^\al) {ai0 


i\Vo) 


= / dX e~ iX ^A x . 

J rp 


(A8) 


Therefore, the constrained amplitude probability A[f\ can be 
written as the sum of the path amplitudes Ay, which are gen¬ 
erated by the effective Hamiltonian Hs{t) — XP (f)A(f). By 
defining the corresponding unitary operator as 

U x (f) = 7V'i ) ( A 9) 


we have 

= {VT\U x (t)\yr 0 ). (A10) 

tP 9 


We can go further with the interpretation of the con¬ 
straint. The parameter X and / can be thought of as eigen¬ 
values of conjugate operators X and / satisfying the rela¬ 
tion (A |/) = e“A/_ Using the relation (^| = 

(A | t we can write 

m = [ ^L(^,ai«-“^-^i<)... 

J rp 

(All) 

Based on Eq. (All), we can interpret A[f \ is the probability 
amplitude to go from the state |y r 0 ,f) to the state A), 
where A and / are conjugate operator acting on an additional 
Hilbert space. The latter is interpreted as the Hilbert space 
of the detector needed to measure the special observable [17— 
19]. The effective Hamiltonian describing the system and a 
quantum detector dynamics is 


H{t) =H s (t) —A P(t)A(t). 


The approach outlined above applies to any time-dependent 
observable A(f). In this work we take A(f) to be the time- 
dependent Hamiltonian. The power operator considered in 
Ref. [15] would be another meaningful choice. 

To determine the variation of the internal energy, we take 
A(f) = Hs(t) and p (f) = 5(T — t) — 5(f). When using a 
discretized evolution, we assume that 5(4 —f) = 1/Af for 
4 < f < 4 + i = 4 + Af and 0 elsewhere. From Eqs. (Bl) we 
have, in the limit of Af —> 0 [18, 19] 

UA7) K, e -i^dH N ~W nAn) e ~mH N - l -Xp N _A N -l) 

)f/(T)e“ ai? (°). (A13) 

Therefore, the total unitary evolution corresponds to two fast 
couplings with the detector with an central driven evolution of 
the system. From this we immediately arrive to the moment 
generating function Gj discussed in the main text. 


Appendix B: System versus environment measurement. 

In the main text we have discussed how to measure the dis¬ 
sipated heat statistics though the system degrees of freedom. 
Here we show that this statistics is the same as we would ob¬ 
tain by measuring directly the environment [3, 35], 

We make the standard assumption that the system and 
the environment are weakly coupled. This allows us to 
neglect the energy related to system-environment coupling 
Hamiltonian. We consider Hs to be time-independent as 
in the fundamental interaction block discussed in the main 
text. The total Hamiltonian reads H = Hs + He + Use ■ 
If we measure the degrees of freedom of the environment, 
we obtain a G x that has the form of G x with U x pJ) = 
exp (iXHe/'X) exp (—iHT) exp (— IXHe /2) [35] . 

In the weak coupling limit, exp (iXHe/2) = 
exp[iX(H — Hs — Hse)/2] « exp [iX(H — Hs)/2] and 

[H . Hs] « 0. Therefore, 

fj x (j) ~ e -aH s /2 e aH/2 e -iHT e -af)/2 e aa s /2 

= e -i*-H s /2 e -iHT e iXH s /2 = jj_^ (j)_ (Bl) 

From Eq. (B 1) it follows that G x = G x and the statistics gen¬ 
erated by measuring He is equal to the one obtained by mea¬ 
suring Hs with opposite sign. The opposite sign in the expo¬ 
nents with respect to Eq. (Bl) takes into account the fact that 
an emission (absorption) by the environment, i.e., decreasing 
(increasing) of the environment energy, corresponds to an ab¬ 
sorption (emission) process of the system, i.e., increasing (de¬ 
creasing) of the system energy. 


Appendix C: First moment of the work done on a quantum 
system 


We first restrict our attention to the dissipated heat and cal¬ 
culate the first moment of the moment generating function G x 


— Q = —i 


t 


■ dGx -Tr \ dUx n f/ + l T! n dU ~ X 

~dX A=o - Tts+e VdX poU -^ UxpQ ^T 


x=o 

(Cl) 


(A12) 
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where U\ = n^ =() II^, U k x = e ~'^fis/ 2 e ~ ,At ^ k e'^s/ 2 , and we 
follow the convention that the heat flowing into the system is 
given a positive sign [3], 

The first term in Eq. (Cl) reads 


■ dUx 
‘ dX 


x=o 


H^U - U n (H$-H^Un-i-Uo +... 


-U N ...U i (H l s -H° S )U 0 - UH° S 


(C2) 


where we used the compact notation U k = e ' A ' /// . In an anal¬ 
ogous way, the second term reads 


dX 


X=o 


1 ■ 
2 ; 




AT-1 


+ u^-K_ 




+ ... 

(C3) 


Putting everything together in Eq. (Cl) and using the cyclic 
property of the trace, it is possible to simplify some of the 
evolution operators U k . After defining p k = U k ...UoPoUQ ...U x 
and po = UqPqUq, we have 


— Q = Tfs \-i; 


-Hs (Pn-i -Pat- 2)--^s(Po-Po) 


• (C4) 


The trace over the system and environment can be 
separated by observing that Trs + # H^ (p N — p^-i) 


Try 




(pat —Pat-i) =Tr s H^(p s ^ N — Ps,n-i) 


. There- 


= Tr<q 

fore, the dissipated heat written in terms of the system degrees 
of freedom reads 


Q = Tr s 


^kHsiPs.k- Ps,k- 1 ) 


= £*Tr s H s Ap sk = Y-kQk 
(C5) 

where Q k is the dissipated heat in the time interval 4_i <t< 


tk- 

Let us now introduce a coupling between the system and 
the detector at the beginning and at the end of the evolution. 
The unitary operator then reads 


U x = e iXA "/ 2 n^U k x e~ iXff s /2 . (C6) 

Accordingly, the calculation in Eq. (C2) is modified as 


if the evolution is unitary, i.e., H k = H k , then the Q k vanish 
and no heat is dissipated. In fact, we have ps. k = U k Ps,k -1 U k , 
[Hs,k,U k } =0, and 


Qk = Tr s 
= Tr, 


H sPs,k~HsPs,k -1 


UkH k s Ps.k-\Ul-H k s p s ^\ 


= 0 . 


(C9) 


The interpretation of the Q k as the dissipative contribution 
to the dynamics is strengthened by the following observation. 
The dissipated heat depends on the variation of the density 
operator A p$ x , which, in turn, can be due to both unitary and 
non-unitary dynamics. However, the unitary contribution to 
the change of p, i.e., the one given by [H$,ps], vanishes iden¬ 
tically when we calculate H^(ps. k — Ps.k - 1 ) ■ Thus, the Q k are 
related solely to the dissipative dynamics. This result is anal¬ 
ogous to the one obtained in Ref. [15]. 

The expressions for Q and W can be written in another 
meaningful way as follows. Instead of grouping A p k , we can 
keep the terms A Hg = H k — fl k 1 as written in Eq. (C2). Then 
the dissipated heat in Eq. (C4) reads 


G = Tr 5 


ps,N + HsPsfi + T. k AHsPs,k 


(CIO) 


As the contributions of the initial and final measurements are 
the same, we have that 


W = Tr s 


EfcA ti^psk 


(Cll) 


This confirms the interpretation discussed in the main text that 
the work can be seen as the instantaneous energy injected in 
the system due to the variation of the Hamiltonian in time. 


Appendix D: Work, heat and entropy in quantum system and 
dynamics 

The Von Neumann entropy in a quantum system is defined 
as S = — 1 Tr[p i slogp,s]. Writing it in the basis {| /)} in which 
Ps is diagonal, we obtain S = — ]T, p s „- log p s „. If we take the 
time derivative of the entropy, we have 

nr = -nPs^EPs.u+PsM) = ^P.v.JogP/, (Dl) 


.dUx(t) 


dX 




0 


A=o 


(W) 


x=o 


(C7) 

which differs from Eq. (C2) by the addition of the term 
HgU - UH° S . 

We find that the average work W is 


W = Trs 


Hs Ps.N - Ps ,0 ~ £ k Hs ( Ps,k - Ps.k- 1 ) 


(C8) 


which differs from Eq. (C5) by the variation of the internal 
energy AH = Try H§ p$.N — H^Psfi ■ Therefore, we have ob¬ 
tained the usual result W = AH — Q [3], 

The heat contributions Q k in (C5) are related to the variation 
of the system density operator A p k during infinitesimal evolu¬ 
tions generated by constant Hamiltonians. We can check that 


where the last equation comes from the fact that JT, p s ii = 0 
due to trace conservation. 

If the decoherence time is smaller than all the system time- 
scale, the system is always in the (time-dependent) thermal- 
ized state. Under this hypothesis the above equation should 
be rewritten explicitly with the time-dependence ps(f) = 
exp [— (iHs(t)]/Zs(t) where /3 = 1 /{ksT) is the inverse tem¬ 
perature of the environment and Zg(t) is the partition func¬ 
tion. In addition, we have log[py(/)] = —Hs{t) — log [Zs{t)\. 
Again, this must be intended in terms of component in the ba¬ 
sis in which ps and //«, are diagonal. If £,■ is the energy of 
the state |i), we have ps = Ps.ii\i){i\ = E,e _ ^ Ei /Zy|i)(tj and 
log ps.u = —) 3e, — log[Zj(f)]. With Eq. (Dl), we can write 

Ll j t = £[/ 3 p« ( e« + \ogZs(t)psM} = Y,Ps.uEi. (D 2 ) 




































Keeping in mind that H$ = Ylk £ k\k){k\, we have that 
L/[ps,//£/] = Tr [^rf-Hs] and we can rewrite the above equation 


as 


Tr 




(D3) 


The variation of the entropy can be related, as in the classical 
case, to the dissipated heat. With this identification we find 


Qk = Tr, 



= k B TASk- 
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